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Abstract. With the aim of deriving the binary frequency 
among Ba and S stars, 56 new spectroscopic orbits (46 and 
10, respectively) have been derived for these chemically- 
peculiar red giants monitored with the CORAVEL spec- 
trometers. These orbits are presented in this paper (38 
orbits) and in a companion paper (Udry et al. 1998, Pa- 
per II; 18 orbits). The results for 12 additional long-period 
binary stars (6 and 6, respectively), for which only min- 
imum periods (generally exceeding 10 y) can be derived, 
are also presented here (10) and in Paper II (2). The global 
analysis of this material, with a few supplementary orbits 
from the literature, is presented in Jorissen et al. (1998). 

For the subsample of Mira S, SC and (Tc-poor) C 
stars showing intrinsic radial- velocity variations due to at- 
mospheric phenomena, orbital solutions (when available) 
have been retained if the velocity and photometric pe- 
riods are different (3 stars). However, it is emphasized 
that these orbit determinations are still tentative. Three 
stars have been found with radial- velocity variations syn- 
chronous with the light variations. Pseudo-orbital solu- 
tions have been derived for those stars. In the case of RZ 
Peg, a line-doubling phenomenon is observed near maxi- 
mum light, and probably reflects the shock wave propa- 
gating through the photosphere. 

Key words: Stars: late-type - Stars: barium - Stars: S - 
Stars: Mira - Binaries: spectroscopic 



1. Introduction 

Several samples of chemically-peculiar red giants (PRG) 
were monitored in radial velocity with the aim of deriv- 
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* Based on observations obtained at the Haute-Provence Ob- 
servatory (France) and at the European Southern Observatory 
(ESO, La Silla, Chile) 



ing their binary frequency. The spectroscopic-binary or- 
bits presented here and in a companion paper (Udry et al. 
1998, Paper II) belong to a systematic survey of barium 
and Tc-poor S stars undertaken to gain insight into the 
formation process of these stars. 

Barium stars and their Pop. II counterparts, the 
CH stars, are not evolved enough to have produced on 
their own the observed overabundances of elements heav- 
ier than Fe (e.g. Lambert 1985), usually associated with 
the nucleosynthesis occurring during He-burning thermal 
pulses on the asymptotic giant branch (AGB). The bi- 
nary nature of these stars accounts for the observed chem- 
ical peculiarities through mass transfer across the binary 
system (McClure et al. 1980; McClure 1983). The exact 
mass-transfer mode is still a matter of debate but is turn- 
ing towards a solution involving Roche-lobe overflow or 
wind accretion, depending on the orbital period (Han et 
al. 1995). 

Contrarily to Tc-rich S stars which are able to synthe- 
size heavy elements on their own, Tc-poor S stars are be- 
lieved to be the cool descendants of barium stars and thus 
also owe their chemical peculiarities to binarity (Jorissen 
& Mayor 1992; Ake 1997; Van Eck et al. 1997). 

The 46 new (43) or updated (3) spectroscopic orbits of 
barium stars (+ 6 Ba stars with P m ; n determinations) and 
the 10 new orbits of S stars (+ 6 binary S stars with P m i n 
determinations) reported in this paper and in Paper II, 
combined with a few other newly available orbits (see ref- 
erences in Jorissen et al. 1998), provide an unequalled 
sample to constrain the binary evolution channels rele- 
vant for barium stars, as identified by Han et al. (1995) 
using only 17 barium-star orbits. It also allows us to re- 
discuss the evolutionary link between barium and Tc-poor 
S stars. This is done in a companion paper (Jorissen et al. 
1998). 

In order not to bias the S-star sample towards low lu- 
minosities, a small sample of S Miras, SC and Tc-poor 
C stars were also monitored. However, the radial-velocity 
jitter associated with atmospheric motions (envelope pul- 
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sation, large convective cells, shocks) induces non-orbital 
radial-velocity variations that limit our ability to detect 
binaries and, a fortiori, to derive their orbital elements, 
at a given level of measurement precision. Illustrative ex- 
amples will be discussed at the end of the paper. 

This paper is organized as follows. Section 2 describes 
the stellar samples. The observations are briefly discussed 
and a short statistical overview of the survey is given in 
Sect. 3. Section 4 provides the orbital elements and the 
radial- velocity curves. Some stars are also commented in- 
dividually. Mira S, SC and C stars with radial-velocity 
variations associated or not with an orbital motion are 
discussed in Sect. 5. 

2. CORAVEL star samples and binary content 

The star samples considered in this paper were observed 
with the CORAVEL spectro-velocimeters (Baranne et al. 
1979) installed on the 1-m Swiss telescope at the Haute- 
Provcncc Observatory (France) and on the 1.54-m Danish 
telescope at the European Southern Observatory (La Silla, 
Chile). The criteria defining these various samples, their 
monitoring history and their binary content are described 
in the following subsections. For a detailed discussion of 
the binary frequency in the corresponding PRG family as 
a whole, we refer to the analysis paper (Jorissen et al. 
1998). 

2.1. Strong barium stars 

This sample includes the 28 barium stars with a strong 
anomaly (Ba4 or Ba5 on the scale defined by Warner 
1965) from the list of Lii et al. (1983), not already moni- 
tored by McClure at the Dominion Astrophysical Observa- 
tory (DAO), plus the special star HD 46407 (see Sect. 2.4). 
These stars are presented in Table 1. The CORAVEL mon- 
itoring of that sample started in 1984 (preliminary re- 
sults were presented in Jorissen & Mayor 1988). All the 
stars show radial-velocity variations due to binary mo- 
tion except HD 19014 for which the situation is still un- 
clear. Among the new spectroscopic orbits of strong bar- 
ium stars, 18 are presented in Sect. 4. The others can be 
found in Paper II (8), in Griffin et al. (1996; 2) and in 
Jorissen et al. (1995; 1). The corresponding references are 
indicated in the last column of the table. 



Table 1. The sample of barium stars with strong anomalies 
(Ba4 or Ba5 on the scale defined by Warner 1965) monitored 
with CORAVEL. The spectral type, visual magnitude and 
B — V index are from Lii et al. (1983) or, if not available, 
from Lii (1991). The last column refers to a note at the end 
of the table providing the reference where the orbital elements 
are to be found 



HD/DM/.. 


Sp. 


Ba 


V 


B — V 


orbit 


5424 


Kl 


4 


9.48 


1.14 


1 


19014 


K4 


5 


8.18 


1.69 





20394 


K0 


4 


8.72 


1.09 


6 


24035 


K4 


4 


8.48 


1.24 


1 


36598 


K2 


4 


8.02 


1.31 


1 


42537 


K4 


5 


8.93 


1.85 


1 


43389 


K2 


5 


8.32 


1.48 


2 


44896 


K3 


5 


7.26 


1.58 


2 


46407° 


K0 


3 


6.28 


1.12 


2 


50082 


K0 


4 


7.42 


1.02 


2 


60197 


K3 


5 


7.74 


1.69 


1 


84678 


K2 


4 


9.00 


1.46 


1 


88562 


K2 


4 


8.54 


1.42 


1 


92626 


K0 


5 


7.10 


1.34 


2 


100503 


K3 


5 


8.74 


1.69 


1 


107541 


K0 


4 


9.39 


1.06 


2 


120620 


K0 


4 


9.62 


1.10 


1 


121447 c 


K7 


5 


7.85 


1.75 


7 


123949 


K6 


4 


8.73 


1.37 


1 


154430 


K2 


4 


8.74 


1.56 


1 


196445 


K2 


4 


9.08 


1.42 


1 


201657 


Kl 


4 


8.04 


1.22 


2 


201824 


K0 


4 


8.91 


1.08 


6 


211594 


K0 


4 


8.11 


1.16 


2 


211954 


K2 


5 


10.21 


1.33 


1 


+38° 118 


K2 


5 


8.88 


1.50 


1 


-42° 2048 


K2 


4 


9.30 


1.41 


1 


-64° 4333 d 


K0 


4 


9.60 


1.20 


1 


Lii 163 


G5 


5 


10.90 


1.18 


1 



a) DAO+CORAVEL measurements 

c) Two classifications: Ba5 or SO. Appears in both tables 

d) CpD 

References to orbital elements: 0: Jorissen et al. (1998); 

1: this paper; 2: Udry et al. (1998, Paper II); 6: Griffin et al. 

(1996); 7: Jorissen et al. (1995) 



2.2. Mild barium stars 

A random selection of 33 stars with a mild barium 
anomaly (Ba<l, Bal, Ba2 on the Warner scale, from the 
list of Lii et al. 1983) has been monitored since 1988 for 
comparison. A few mild barium stars observed by the Mar- 
seille team (see Paper II) for a different purpose were in- 
cluded in our initial sample later on. The global sample, 
also including 3 DAO stars (see next section), is presented 
in Table 2. 



Among these mild barium stars, 27 stars are definitely 
binaries (3 with fixed parameters, and 6 with only a lower 
limit available on the orbital period) , 3 are suspected bina- 
ries, 4 show no evidence of binary motion and 2 (HD 65699 
and HD 206778) are supcrgiants misclassified as mild bar- 
ium stars (Smith & Lambert 1987; superscript in Ta- 
ble 2). Results for 14 mild barium stars are presented in 
Sect. 4. The remaining 13 binaries are to be found in Pa- 
per II (12) and in Griffin (1991; 1). 
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Table 2. The sample of barium stars with mild anomalies 
(Ba<l, Bal, Ba2 on the scale defined by Warner 1965) mon- 
itored with CORAVEL. The spectral type, visual magnitude 
and B — V index are from Lii et al. (1983) or, if not available, 
from Lii (1991). Column 6 gives a binarity flag (o: orbit, po: 
preliminary orbit with fixed parameters, mp: minimum period 
only available, sb: suspected binary, c: no radial-velocity vari- 
ation). The last column refers to a note at the end of the table 
providing the reference where the orbital elements are to be 
found 

HP/DM Sp. Ba V B-V bin. o7bit~ 
18182 KO < 1 8.95 1.01 sb 

22589 G5 < 1 9.00 0.72 po 1 

26886 G8 1 7.97 0.97 o 2 

27271 G8 1 7.49 1.00 o 2 

40430 K0 1 8.07 1.01 mp 1 

49841 G8 1 8.56 0.99 o 2 

50843 Kl 1 8.13 1.06 c 

51959 K2 1 8.95 1.06 mp 1 

53199 G8 2 9.08 0.95 po 2 

58121 K0 1 7.93 1.14 o 2 

59852 G9 1 8.67 0.92 po 1 

65699° G5 1 5.07 1.12 
91208 K0 1 8.02 0.95 o 1 

95193 K0 1 8.28 0.99 o 1 

95345 K2 1 4.83 1.16 c 

101079 Kl 1 8.43 0.96 mp 2 
104979 K0 1 4.16 0.99 mp 2 
119185 K0 1 8.91 1.00 c 
130255 K0 1 8.61 1.08 c 
131670° Kl 1 8.00 1.20 o 1 
134698 Kl 1 8.72 1.31 mp 1 
139195 Kl 1 5.26 0.95 o 12 
143899 G8 1 8.30 1.08 o 1 
180622 Kl 1 7.62 1.23 o 2 
183915 K0 2 7.28 1.34 sb 
196673° K2 2 6.98 1.12 o 1 
200063 K3 1 7.31 1.60 o 2 
206778° K2 < 1 2.38 1.52 
210946 Kl 1 8.07 1.08 o 2 
216219 G5 1 7.45 0.74 o 2 
218356 K2 2 4.47 1.35 sb 
223617° K2 2 6.94 1.16 o 2 
288174 K0 1 8.99 1.24 o 1 
-01°3022 Kl 1 9.25 1.09 o 1 
-10°4311 GO 1 9.89 0.63 mp 1 
-14°2678 K0 < 1 9.81 0.99 o 1 

a) DAO+CORAVEL measurements 

b) misclassified as mild Ba star 

References to orbital elements: 0: Jorissen et al. (1998); 
1: this paper; 2: Paper II; 12: Griffin (1991) 



2.3. DAO stars 

For several barium stars monitored by McClure, a few re- 
cent CORAVEL measurements, obtained in the framework 
of other projects, are nevertheless available. These new 



measurements often allow us to improve the DAO orbit, 
since they significantly increase the span of the monitor- 
ing. All these new orbital parameters are used in the analy- 
sis paper (Jorissen et al. 1998). The present paper and Pa- 
per II provide only four of these stars, 3 for which the num- 
ber of CORAVEL measurements is fairly large (HD 46407, 
HD 131670 and HD 223617) and 1 (HD 196673) for which 
the 2 new available observations significantly change the 
period obtained by McClure & Woodsworth (1990). These 
stars are also included in Tables 1 and 2, where they are 
identified by superscript 'a' in column 1. Their orbits are 
presented in Sect. 4 or in Paper II. 



2.4- Two barium stars of special importance 

The K0Ba3 star HD 46407 and K2Ba2 star HD 218356 
(=56 Peg) were included in the CORAVEL samples (Ta- 
bles 1 and 2) because of their unique photometric be- 
haviour among barium stars. HD 46407 exhibits long-term 
photometric variations in phase with the orbital motion 
(Jorissen et al. 1991; Jorissen 1997), whereas 56 Peg is a 
strong X-ray and UV source, indicating that it is an in- 
teracting binary system possibly hosting an accretion disk 
(Schindler et al. 1982; Dominy & Lambert 1983). A com- 
bined DAO / CORAVEL orbit of HD 46407 is presented in 
Paper II whereas no clear orbital solution emerges for 56 
Peg yet. 



2.5. Non-variable S stars 

An initial sample of 9 S stars, whose monitoring started 
in 1984 (with preliminary results presented in Jorissen & 
Mayor 1988, 1992) was extended to 36 stars in 1988. This 
sample contains bright northern S stars from the General 
Catalogue of Galactic S Stars (GCGSS; Stephenson 1984) 
with no variable star designation, neither in the General 
Catalogue of Variable Stars (Kholopov et al. 1985) nor in 
the New Catalogue of Suspected Variable Stars (Kukarkin 
et al. 1982). The criterion of photometric stability has 
been adopted to avoid the confusion introduced by the 
envelope pulsations masking the radial-velocity variations 
due to orbital motion. Such a selection criterion clearly 
introduces a strong bias against intrinsically bright S stars. 

The sample of 36 photometrically non- variable S stars 
is presented in Table 3. Among them, 24 are binaries 1 (6 
have only a lower limit available on the orbital period), 
10 show no evidence for binary motion, and 2 (HD 262427 
and BD +22°4385) are likely misclassified (Jorissen et al. 
1998). The new spectroscopic binaries (18) are described 
in Sect. 4 (16) and in Carquillat et al. (1998; 2). 



1 including HD 121447 already presented in the strong bar- 
ium sample 
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K0 
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8.95 


1.01 


sb 


22589 


G5 


< 1 


9.00 


0.72 


po 


26886 


G8 


1 


7.97 


0.97 
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27271 
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7.49 


1.00 


o 


40430 


K0 
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8.07 


1.01 


mp 


49841 


G8 


1 


8.56 


0.99 


o 


50843 


Kl 


1 


8.13 


1.06 


c 


51959 


K2 


1 


8.95 


1.06 


mp 


53199 


G8 


2 


9.08 


0.95 


po 


58121 


K0 


1 


7.93 


1.14 


o 


59852 


G9 


1 


8.67 


0.92 


po 


65699° 


G5 


1 


5.07 


1.12 




91208 


K0 


1 


8.02 


0.95 


o 


95193 


K0 


1 


8.28 


0.99 


o 


95345 


K2 


1 


4.83 


1.16 


c 


101079 


Kl 


1 


8.43 


0.96 


mp 


104979 


K0 


1 


4.16 


0.99 


mp 


119185 


K0 


1 


8.91 


1.00 


c 


130255 


K0 


1 


8.61 


1.08 


c 


131670° 


Kl 


1 


8.00 


1.20 


o 


134698 


Kl 


1 


8.72 


1.31 


mp 


139195 


Kl 


1 


5.26 


0.95 


o 


143899 


G8 


1 


8.30 


1.08 


o 


180622 


Kl 


1 


7.62 


1.23 





183915 


K0 


2 


7.28 


1.34 


sb 


196673° 


K2 


2 


6.98 


1.12 





200063 


K3 


1 


7.31 


1.60 


o 


206778° 


K2 


< 1 


2.38 


1.52 




210946 


Kl 


1 


8.07 


1.08 





216219 


G5 


1 


7.45 


0.74 





218356 


K2 


2 


4.47 


1.35 


sb 


223617° 


K2 


2 


6.94 


1.16 


o 


288174 


K0 


1 


8.99 


1.24 


o 


-01°3022 


Kl 


1 


9.25 


1.09 


o 


-10°4311 


GO 


1 


9.89 


0.63 


mp 


-14° 2678 


K0 


< 1 


9.81 


0.99 


o 
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Table 3. The sample of photometrically non-variable S stars 
monitored with CORAVEL. The spectral type is from the 
GCGSS. The V magnitude and B — V index for S stars are from 
various sources, as listed in The General Catalogue of Photo- 
metric Data (GCPD; Mermilliod et al. 1997). In case no data 
are listed in the GCPD, the V magnitude from the GCGSS is 
listed. Column 6 gives a binarity flag (o: orbit, po: preliminary 
orbit with fixed parameters, mp: minimum period only avail- 
able, j: jitter, c: no radial- velocity variation). The last column 
refers to a note at the end of the table providing the reference 
where the orbital elements are to be found 



HD/DM 


GCGSS 


Sp. 


V 


B—V 


bin. 


orbit 


7351 


26 


S3/2 


6.43 


1.70 


o 


3 


30959 


114 


bo/ 1 


4. (4 


1 OA 

1.84 


mp 


i 
1 


35155 


133 


S4,l 


6.77 


1.80 


o 


0,10 


49368 


260 


S3/2 


7.65 


1.78 


o 


1 


61913 


382 


M3S 


5.56 


1.64 


j 





63733 


/iii 
411 


O A /O 

S4/3 


t n a 

7.94 


1 7/1 

1.74 


o 


1 


95875 


720 


S3, 3 


8.74 


1.11 


o 


1 


1 1 1 A A 7C 

121447 




OA 
SO 


7.85 


1.75 





1 


170970 


1053 


S3/1 


7.60 


2.00 





1 


184185 


1140 


S3*4 


9.20 


2.00 


mp 


1 


lo9ool 


1 1 *7Q 
11(8 


b4 Z 


o. 10 


o nn 
z.UU 


j 


n 
U 


191226 


1192 


M1-3S 


7.34 


1.32 





3 


i ni con 


i id/I 
1194 


b 


/. tu 


o nn 
Z.UU 





i 
1 


1 AO A ACt 

iy244o 


1198 


So/1 


9.80 


2.00 


j 





2100 * 2 


1 oi r 

lolo 


Q/1 /I 


b.oo 


i Qn 
l.oU 


j 


n 
U 


21o0o4 


1 QOO 

lozz 




k no. 


1 A 7 

1.4/ 


mp 


1 


24oolo 


156 


s 


9.60 


2.00 


o 


1 


262427 


247 




9.90 


2.00 







2ooooo 


Zoo 


CO /o 
Do/ Z 


o a n 
y.4U 


i.yi 


mp 


1 


332077 


1201 


S3,l 


9.00 


1.96 


o 


10 


343486 


1092 


S6,3 




2.00 


o 


1 


+04°4354 


1193 


S4*3 


9.50 


2.00 


j 





+15° 1200 


219 


S4/2 


9.39 


2.01 


j 





+20°4267 


1158 


Swk 


10.80 


2.55 


j 





+21° 255 


45 


S3/1 


8.60 


1.30 





1 


+22° 700 


96 


S6,l 


10.50 


2.00 


o 


10 


+22°4385 6 


1271 




10.10 


2.00 







+23°3093 


981 


S3/3 


9.90 


2.00 





10 


+23°3992 


1209 


S3,3 


10.50 


2.00 





1 


+24° 620 


87 


S3/3 


8.88 


2.06 





10 


+28°4592 


1334 


S2/3: 


9.50 


1.40 


o 


1 


+31°4391 


1267 


S2/4 


9.30 


2.00 


mp 


1 


+79° 156 


106 


S4/2 


10.10 


1.70 


mp 


1 


-04°2121 


416 


S5/2 


9.20 


2.00 


j 





-10° 1334 


176 


Sr 


9.00 


2.00 


j 





-21°2601 


554 


S3*3 


9.20 


2.00 


j 






Table 4. The sample of Mira S stars, SC stars and C stars with 
no Tc lines, monitored with CORAVEL. The spectral type is 
from the GCGSS, and B — V from the GCPD 



HD/DM 


GCGSS 


Var 


Sp. 


B — V 


1967 


9 


R And 


S5-7/4-5e 


1.97 


4350 


12 


U Cas 


S5/3e 


2.00 


14028 


49 


W And 


S7/le 


1.70 


29147 


103 


T Cam 


S6/5e 


2.30 


53791 


307 


R Gem 


S5/5 


2.10 


70276 


494 


V Cnc 


S3/6e 


2.20 


110813 


803 


S UMa 


S3/6e 


2.10 


117287 


— 


R Hya 


M6e-M9eS 


1.60 


185456 


1150 


R Cyg 


S6/6e 


2.00 


187796 


1165 


X Cyg 


S7/1.5e 


1.80 


190629 


1188 


AA Cyg 


S6/3 


1.80 


195763 


1226 


Z Del 


S4/2e 


2.00 


211610 


1292 


X Aqr 


S6,3e: 


2.00 


286340 


117 


GP Ori 


SC7/8 


2.90 


44544 


212 


FU Mon 


S7/7 


3.00 


-04°1617 


244 


V372 Mon 


SC7/7 


2.50 


-08° 1900 


344 




S4/6 


1.90 


54300 




R CMi 


CS 


2.50 


198164 




CY Cyg 


CS 


2.00 


209890 




RZ Peg 


CS 


3.90 


46687 




UU Aur 


C no Tc 


3.00 


76221 




X Cnc 


C no Tc 


3.36 


108105 




SS Vir 


C no Tc 


3.00 



low- luminosity S stars (see Sect. 2.5). However, the enve- 
lope pulsations of Mira stars seriously hamper that search 
by causing a substantial radial- velocity jitter. This will be 
discussed in Sect. 5. 

A sample of 7 SC and CS stars suffering the same 
problem has been monitored as well with CORAVEL, along 
with the 3 carbon stars lacking Tc from the list of Little 
et al. (1987). These stars are also listed in Table 4. Radial- 
velocity variations are observed for all these stars. 

In case the photometric and radial- velocity periods are 
different, a tentative orbital solution is proposed. The re- 
sults are presented in Sect. 5. 

3. Observations and overview of the survey 

Since 1983, more than 2500 measurements have been gath- 
ered for the stars of the samples described above. The ra- 
dial velocities arc obtained by usual cross-correlation (cc) 
and gaussian-fitting techniques (see Baranne et al. 1979 
or Duquennoy et al. 1991 for details). An illustrative ex- 
ample of a profile obtained for HD 24035 (cc-dip with the 
fitted gaussian) is displayed in Fig. la. 

The still unpublished individual measurements will 
be available at the Centre de Donnees Stellaires 
(CDS) in Strasbourg or on our dedicated web page 
(obswww.unige.ch/^udry/cine/barium/barium.html) 
and thus will not appear here. 



h) star misclassified as S? 

c) Two classifications: Ba5 or SO. Appears in both tables 
References to orbital elements: 0: Jorissen et al. (1998); 1: this 
paper; 3: Carquillat et al. (1998); 10: Jorissen & Mayor (1992) 



2. 6. Mira S stars, SC stars and Tc-poor carbon stars 

A sample of 13 Mira S stars (Table 4) has also been mon- 
itored in order not to restrict the search for binaries to 
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Fig. 1. Statistical overview of the combined samples of bar- 
ium and S stars: a) typical cross-correlation profile (here an ob- 
servation of HD 24035); b) distribution of individual measure- 
ment errors Si (in kms -1 ); c) distribution of the number of 
CORAVEL measurements per star N; d) distribution of spans 
in days (AT) between the first and the last measurements, for 
every star 

The mean precision of the measurements is about 
0.3 kms -1 (more than 76% of the measurements be- 
tween 0.25 and 0.4 kms -1 ) as shown in Fig. lb display- 
ing the distribution of individual errors £j. S stars with 
non-orbital radial- velocity variations usually present non- 
gaussian cc-dips (see Sect. 5.1) leading to non-realistic er- 
ror estimates that slightly swell the upper tail of the dis- 
tribution. Stars with broad cc-dip widths act in the same 
way. 

The distributions of the number of observations and 
of the time span of the observations for every star are 
given in Figs, lc and Id, respectively. They illustrate the 
large observational effort devoted to these programmes. 
The median value for the number of observations per star 
is N = 19 which permits a good estimate of the orbital 
parameters. The typical span (AT = 3398 d) is about 10 
years ensuring a good completeness of the orbital periods 
up to this value. 

Further information on the CORAVEL observation and 
reduction techniques can be found in Duquennoy et al. 
(1991). 

4. Radial-velocity curves and orbital parameters 

The radial-velocity monitoring of a binary star may lead 
to different qualitative results depending on the parame- 



ters of the orbital solution (period P, eccentricity e, ampli- 
tude K), the number of measurements or the observation 
sampling. The strategy for the presentation of the results 
will thus depend on these different cases. Situations pos- 
sibly encountered are the following: 

1) A stable solution is found and the orbital parame- 
ters can be derived (flag 'o' in Tables 2 and 3). This is the 
case for most of the orbits presented in this paper. A table 
is provided with the orbital parameters and their uncer- 
tainties and a few additional interesting related quantities 
as the number of measurements (N) used to derive the or- 
bital solution or the residue (O — C) around this solution. 
The phase-folded orbital solution is then also displayed. 
Badly constrained parameters are readily identifiable by 
their large uncertainties. 

2) The orbital solution exists but is not fully con- 
strained. One or several parameters (usually P or e) have 
to be fixed (e.g. for an uncompletely-covered orbit or 
when the periastron passage in an eccentric orbit has been 
missed; flag 'po' in Tables 2 and 3). In such a case, the 
adequate orbital parameter is fixed to a probable value as 
given by the orbital solution with minimum residuals. The 
obtained minimized orbital elements are given in the same 
table as for case 1) but with no uncertainties on the fixed 
parameter (s). A diagram with the velocity measurements 
folded in phase is also provided. 

3) The star has a clearly variable radial velocity but 
the period is insufficiently covered (usually just a drift is 
observed) to derive a preliminary orbit, even with fixed 
parameters (flag 'mp' in Tables 2 and 3). No solution is 
found. The minimum period (P > P m i n ) is indicated in 
the table of orbital elements, along with the number of 
measurements and the time span of the observations. The 
figure only displays radial velocities as a function of Julian 
date. In addition to the individual measurements available 
at the CDS (Sect. 3), the analysis paper (Jorissen et al. 
1998) summarizes the interesting averaged quantities (ra- 
dial velocity and corresponding uncertainty, etc.) for these 
stars. 

4) No orbital solution can be found because of other 
sources of radial- velocity variations (jitter, flag 'j' in Ta- 
ble 3) masking a possible orbital motion. Average quanti- 
ties are given in Jorissen et al (1998). 

In the following subsections the results for the vari- 
ous star samples described in Sect. 2 will be presented and 
discussed in turn. 

4-1- Strong barium stars 

Definitive or preliminary orbits have been obtained for all 
strong barium stars except HD 19014. The orbital param- 
eters are given in Table 5. For two uncompletely-covered 
orbits (HD 123949 and HD 211954), the periods have been 
fixed to minimized values (see item 2 above). In 8 cases, 
even though the Lucy-Sweeney test was compatible with 
a circular orbit at a 5% confidence level (Lucy & Sweeney 




Fig. 2. Phase-folded radial-velocity curves for the strong barium stars. Minimized periods (see item 2 in Sect. 4) were fixed for 
HD f 23949 and HD 211954 because of the non-complete coverage of the orbits. BD+38°118 is a triple hierarchical system 
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Table 5. Orbital elements for strong barium stars. No uncertainties are given for fixed parameters. The symbol > is used for 
uncertainties exceeding the parameter values in case of badly constrained orbits. N is the number of measurements used to 
derive the orbital solution and O — C the residue around this solution. AT is the span of the observations 
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Fig. 3. Phase-folded radial- velocity curves for mild barium stars with orbital solutions. Open circles are for DAO measurements. 
Two new CORAVEL measurements allow us to propose a new minimum-period estimate for HD 196673. Long-period stars 
without orbital solutions have their radial velocities displayed as a function of Julian dates 
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Fig. 3. (continued) 



1971), the slightly 'eccentric' orbit has been listed, be- 
cause in the case of barium stars, there is no physical argu- 
ment to prefer the circular orbit (see Jorissen et al. 1998). 
The corresponding radial- velocity measurements folded in 
phase are displayed in Fig. 2. 



BD +38° 118 is a triple hierarchical system. The long- 
period orbit (noted ab+c) describes the motion of the c 
component, relatively to the center of mass of the close 
(short-period) system (noted a+b). The orbits are ob- 
tained iteratively by correcting the short-period orbital 
motion from the long-period perturbation. 



4-2. Mild barium stars 

Orbital elements were derived for 10 among the 14 mild 
barium stars presented in this section (Table 6), the others 
only allowing minimum period estimates. The correspond- 
ing phase-folded curves are displayed in Fig. 3, along with 
the temporal radial- velocity variations of the remaining 4 
stars with no orbital solution. 

The star HD 196673 deserves a special note. It belongs 
to the DAO sample with a spectroscopic orbit published 
by McClure & Woodsworth (1990). However, as shown 
by 2 new CORAVEL measurements, the inferred orbital 
period (P = 4000 d) was too short. Based on all the 
DAO+CORAVEL measurements, we propose a new period 




Fig. 4. Phase-folded radial- velocity curves of the photometrically non-variable S stars with orbital solutions. Long-period stars 
without orbital solutions have their radial velocities displayed as a function of Julian dates 
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Fig. 4. (continued) 
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Table 6. Orbital elements for mild barium stars. No uncertainties are given for fixed parameters. N is the number of measure- 
ments used to derive the orbital solution and O — C the residue around this solution. AT is the span of the observations 



Id 


P 


T [HJD 


e 


7 




K 


f(m) 


asini 


N 


O-C 


AT 


HD/DM 


[days] 


-2400000] 




r i — li 

[kms ] 


[deg] 


r i — li 

[kms ] 




[Gm] 




r i — li 

[kms ] 


[days] 


22589 


5721.181 


48674.73 


0.240 


-27.97 


341.09 


1.97 


4.180e-03 


150.74 


19 


0.217 


3641 




453.986 


133.44 


0.175 


0.62 


17.38 


0.28 


2.471e-03 


63.26 








40430 


P > 3700 
















15 




3720 


51959 


P > 3700 
















21 




3720 


59852 


3463.906 


46841.03 


0.152 


0.07 


88.83 


1.85 


2.190e-03 


86.98 


19 


0.273 


3643 




53.776 


187.14 


0.057 


0.07 


20.68 


0.12 


4.376e-04 


5.93 








91208 


1754.002 


45628.36 


0.171 


0.19 


79.25 


5.05 


2.239e-02 


119.93 


24 


0.387 


3036 




13.321 


54.16 


0.022 


0.08 


10.11 


0.11 


1.555e-03 


2.91 








95193 


1653.728 


46083.62 


0.135 


-7.32 


283.23 


5.37 


2.585e-02 


120.97 


18 


0.264 


2596 




8.950 


38.25 


0.017 


0.07 


7.60 


0.10 


1.446e-03 


2.34 








131670 


2929.730 


46405.11 


0.162 


-25.13 


49.88 


5.16 


4.005e-02 


204.94 


55 


0.363 


5865 




12.265 


44.96 


0.014 


0.05 


5.44 


0.08 


1.830e-03 


3.22 








134698 


P > 3600 
















22 




3630 


143899 


1461.608 


46243.43 


0.194 


-29.83 


277.66 


4.95 


1.740e-02 


97.64 


26 


0.271 


2649 




6.919 


19.76 


0.025 


0.06 


4.57 


0.10 


1.065e-03 


2.04 








196673 


6500.0 


43953.72 


0.645 


-24.61 


111.06 


3.49 


1.274e-02 


237.99 


51 


0.468 


5971 




— 


26.79 


0.028 


0.10 


4.48 


0.17 


2.175e-03 


13.54 








288174 


1824.286 


47157.62 


0.194 


34.67 


48.33 


4.55 


1.684e-02 


111.96 


14 


0.147 


3296 




7.057 


17.73 


0.015 


0.04 


3.95 


0.07 


7.431e-04 


1.70 








-01°3022 


3252.530 


46817.14 


0.285 


-35.38 


212.09 


3.76 


1.586e-02 


161.34 


26 


0.253 


3667 




31.420 


34.35 


0.018 


0.06 


4.31 


0.09 


1.185e-03 


4.28 








-10° 4311 


P > 3400 
















33 




3392 


-14° 2678 


3470.519 


48828.06 


0.217 


4.93 


282.61 


4.07 


2.252e-02 


189.38 


15 


0.390 


3295 




107.690 


107.22 


0.040 


0.12 


12.16 


0.14 


2.501e-03 


8.94 









(fixed to P = 6500 d), which is a lower bound to the actual 
period. The new preliminary orbital elements obtained by 
fixing the period to the above value are listed in Table 6. 



4-3. Photometrically non-variable S stars 

Among the 16 binary S stars with no photometric varia- 
tions presented in this subsection, an orbital solution has 
been derived for 10 of them whereas 6 have only minimum- 
period estimates. The results are given in Table 7. The 
corresponding phase diagrams are displayed in Fig. 4. 

Note that BD +21°255 is a visual binary, most prob- 
ably of optical nature (Jorissen & Mayor 1992). Table 7 
and Fig. 4 provide the orbit of the S star (BD +21°255 
= PPM 91178 = SAO 75009 = HIC 8876), whereas the 
visual K-type companion (BD +21°255p = PPM 91177 = 
SAO 75008) is also a spectroscopic binary whose orbit is 
given in Jorissen & Mayor (1992). 

The jitter level of the non-binary stars of the sample 
is of the order of 1-2 kms -1 (Jorissen et al. 1998). Its 
hampering influence on the detection of binarity is thus 
limited to low-amplitude orbital motions. 



5. Mira S stars, SC and C stars: intrinsic radial- 
velocity variations versus orbital motion 

5.1. Intrinsic radial-velocity variations of Mira variables 

The large velocity jitter observed in Mira S stars (Ta- 
bles 3d, 3e and 4 in Jorissen et al. 1998) is related to 
their complex and variable cc-dips (Barbier et al. 1988). 
In some cases like for the star \ Cyg, the cc-dip has an in- 
verse P-Cygni shape (Fig. 5a), out of which a meaningful 
radial velocity is very difficult to extract. In other cases 
like AA Cyg and R Hya, the cc-dips are featureless, broad 
and very stable, and these stars have the smallest jitter 
in our sample of Mira S stars (Table 3d of Jorissen et 
al. 1998). For stars like R And (=HD 1967) and RZ Peg 
(=HD 209890), two clearly distinct dips are present. In 
RZ Peg for example (Figs. 5c-d), the two dips correspond 
to velocities of —25 and —50 kms -1 , whereas the center 
of mass of the star (as probed by submm observations of 
the circumstellar CO rotational lines) moves with a ve- 
locity of —35 kms -1 (Sahai & Liechti 1995), close to the 
average of the two CORAVEL dips. The two minima ob- 
served in the cc-dip are therefore likely associated with 
upwards- and downwards-moving layers forming a shock 
in the Mira atmosphere (e.g. Fox et al. 1984, Querci 1986, 
Bowcn 1988, Bcsscll et al. 1996). The respective weights 
of the two minima vary with time, and so does the ve- 
locity of the resulting blend to which a gaussian is fitted. 
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Table 7. Orbital elements for the S stars with no (strong) light variations. The symbol > is used for uncertainties exceeding 
the parameter values in case of badly constrained orbits. N is the number of measurements used to derive the orbital solution 
and O — C the residue around this solution. AT is the span of the observations 
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Table 8. Orbital or pseudo-orbital (pulsational) elements for Mira S stars, SC stars and Tc-poor C stars. The most likely cause 
of the radial-velocity variations is given in column 2: 'orb' for orbital motion and 'puis' for intrinsic atmospheric phenomenon. 
No uncertainties are given for fixed parameters. The symbol > is used for uncertainties exceeding the parameter values in case 
of badly constrained orbits. N is the number of measurements used to derive the orbital solution and O — C the residue around 
this solution. AT is the span of the observations 
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As an illustration, 2 profiles at different phases (</W) are 
given in Figs. 5c-d for RZ Peg. The derived mean radial 
velocities are indicated on the figures. 



When the CORAVEL cc-dip of RZ Peg is fitted by a 
single gaussian profile, as indicated on Fig. 5, a pseudo- 
orbital solution is found (Table 8), with a period of 
437.3±3.9d. This radial-velocity period may be identified 
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HD 110813 S UMa orb HD 76221 X Cnc orb? 




4.7xl0 4 4.8xl0 4 4.9xl0 4 5xl0 4 4.7xl0 4 4.8xl0 4 4.9xl0 4 5xl0 4 

Julian Date - 2400000 Julian Date - 2400000 

Fig. 6. Phase- folded radial-velocity curves of Mira, SC and S stars with tentative- or pseudo-orbital solutions. Labels ('orb', 
'orb?' or 'puis') are used to characterize our confidence level in the obtained solution or to indicate intrinsic radial-velocity 
variations. Two stars (BD —08° 1900 and HD 286340) have also their radial velocities displayed as a function of Julian dates 
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Fig. 5. Examples of cross-correlation dips for Mira S stars, 
a-b) Profiles for \ Cyg, including an inverse P Cygni-like pro- 
file; c-d) Two profiles at different phases (</>rv) for the star 
RZ Peg. The time- varying weights of the two minima induce a 
variation of the measured radial velocity 



with the photometric period of 437.8 d derived from the 
AAVSO light curve for the Mira variations (see Sect. 5.2). 
It is therefore clear that, although they mimick an orbital 
motion, the radial-velocity variations of RZ Peg arc in- 
trinsic to the Mira phenomenon. They will be discussed 
in more details in Sect. 5.2. The same coincidence be- 
tween the periods of the light variations and of the ve- 
locity variations is obtained for the CS Mira star R CMi 
(=HD 54300) and for the C semi-regular (SRa) star SS Vir 
(=HD 108105) as well. The radial-velocity variations of 
R CMi can be fitted with a pseudo-orbital solution of pe- 
riod 337.3±1.3d, whereas the GCVS lists P = 337.8 d for 
the associated Mira variations. It is interesting to note 
that the radial-velocity variations of R CMi are strongly 
non-sinusoidal in spite of the almost symmetric light curve 
(the rise time from minimum to maximum light repre- 
sents 48% of the total cycle, according to the GCVS, com- 
pared to 44% for RZ Peg). For SS Vir, the radial- velocity 
and photometric periods are 361.2±3.4d and ~ 364 d, re- 
spectively. The relevant parameters of the pseudo-orbits 
for RZ Peg, R CMi and SS Vir are given in Table 8. The 
radial-velocity periods have been used to fold the points 
in phase. The corresponding diagrams are displayed in 
Fig. 6. 



5.2. Intrinsic radial-velocity variations of Mira S stars: 
the cases of RZ Peg and \ Cyg 

The radial-velocity curve of RZ Peg is well sampled, and 
makes it possible to investigate in more details the prop- 
erties of the double cc-dip phenomenon. 

Figure 7 presents the evolution of the CORAVEL cross- 
correlation profiles as a function of the radial-velocity 
phase (derived from the elements listed in Table 8). The 
evolution of the profiles, from a single cc-dip to a dou- 
ble cc-dip and vice versa, is quite obvious. The radial ve- 
locities corresponding to each of the peaks have been ex- 
tracted by fitting a double-gaussian function to the profile. 
The resulting velocity curve as a function of the phase is 
displayed in Fig. 8. Two conclusions may be drawn from 
this figure: (i) each peak has a rather constant velocity, 
with the center-of-mass velocity corresponding almost to 
the average of the two peaks, and (ii) the double dip oc- 
curs around velocity phases 0rv = 0.7 - 0.9, and this 
behaviour repeats from one cycle to the other. The veloc- 
ity difference between the two peaks is of the order of 20 
- 30 kms" 1 . 

The light curve of RZ Peg has been kindly provided to 
us by the AAVSO (Mattei 1997), and makes it possible to 
convert velocity phases into photometric phases. The pho- 
tometric period appears reasonably stable over the time 
span covered by the radial-velocity observations, with an 
average period of 437.8 d, very close to the 437.3 d period 
of the velocity variations. Adopting a period of 437.3 d, the 
relation <f) p hot = 4>RV + 0.26 is derived from the AAVSO 
maximum (HJD 2 448 301) closest to the zero point of the 
radial-velocity phase (HJD 2 448 413). From this relation, 
it may be concluded that the double cc-dip occurs just af- 
ter maximum light (0 p hot = 0.0 — 0.2). 

A similar phenomenon of line doubling near maximum 
light has been reported in several other Mira stars. In the 
Mira S star v Cyg, Hinkle et al. (1982) reported such a 
line doubling for the rotation-vibration second overtone 
(Av = 3) infrared CO lines and high-excitation first over- 
tone (Av = 2) CO lines. Maehara (1971) also reported the 
doubling of CN and atomic lines near 800.0 nm around 
maximum light in \ Cyg, and Gillet et al. (1985) in S 
Car. This line-doubling phenomenon is the signature of 
a shock front propagating through the line-forming re- 
gion (e.g. Gillet et al. 1985, Bessell et al. 1996). Inci- 
dentally, in x Cyg (Fig. 5), a line-doubling as clear as in 
RZ Peg is not observed for the blue-violet iron lines sam- 
pled by CORAVEL (Baranne et al. 1979). The CORAVEL 
radial velocities of \ Cyg (averaging to +3.2 kms -1 ; Ta- 
ble 3d of Jorissen et al. 1998), derived from the clean cc- 
dips (Fig. 5b), agree with the value obtained by Pierce 
et al. (1979). Hinkle et al. (1982; their Fig. 15) argue that 
these blue- violet photosphcric absorption lines with a con- 
stant radial velocity form in the infalling material for the 
whole pulsation cycle. Nevertheless, the CORAVEL profile 
of x Cyg exhibits some characteristic features that may be 
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Fig. 7. CORAVEL cc-dips of RZ Peg as a function of the 
radial-velocity phase (0 RV )- Phase (HJD 2 448 413.3) is close 
to the epoch of minimum radial- velocity (Fig. 6). The profiles 
have been fitted with a single or double gaussian function, when 
appropriate (thin solid line) 



related to the shock front propagating through the photo- 
sphere. Near visual phases 0.0-0.2, the CORAVEL profile 
of x Cyg has a typical inverse P-Cygni shape (Fig. 5a) 
with the blue-shifted 'emission' component peaking at 
—20 kms -1 . This velocity corresponds to that of the out- 
flowing material, as derived from the second overtone CO 
rotation- vibration lines by Hinkle et al. (1982). A simi- 
lar inverse P-Cygni profile has been reported by Ferlet 
& Gillet (1984) for the Til lines near 1 /im in Mira near 
maximum light. 

It is very likely that the same physical phenomenon 
- namely a shock front propagating through the line- 
forming region - is responsible for the time-dependent 
features observed in the CORAVEL cc-dips of RZ Peg and 
X Cyg. This conclusion is suggested by the fact that in 
both stars, the line-doubling occurs near maximum light, 
and the offset between the two distinct components of 
the cc-dip is of the same order (20 to 30 kms" 1 , corre- 
sponding to the shock discontinuity) . These observed fea- 
tures are well reproduced by the synthetic Fel and CO line 
profiles computed by Bessell et al. (1996) in a dynamical 
Mira atmosphere with a propagating shock. For some rea- 




Fig. 8. Radial velocities of RZ Peg as a function of the pho- 
tometric phase (<j} p hot = 4>rv + 0.26). Open and filled circles 
indicate the velocities of the double-line reductions whereas 
star-like symbols are used for velocities derived from gaus- 
sian one-line fitting. The center-of-mass velocity, as probed by 
submm observations of the CO rotational lines (Sahai & Liechti 
1995), is represented by the horizontal dashed line 



son however, the atmospheric structures of RZ Peg and 
X Cyg must be different, so that only in RZ Peg are the 
blue- violet iron lines sampled by CORAVEL forming in ab- 
sorption on both sides of the shock, resulting in a clean 
double- minima cc-dip. 

5.3. Binaries among Mira variables and SC stars 

Because they are difficult to observe at minimum light, 
Mira S and C stars were generally insufficiently sampled 
to attempt a period search on their radial-velocity curve. 
Well-sampled curves are, however, available for the Mira 
S stars AA Cyg, % Cyg, R Hya, for the C (no Tc) stars 
SS Vir and S UMa, and for all the SC stars. Some of these 
stars have already been discussed in Sects. 5.1 and 5.2. 

No satisfactory periods emerge for AA Cyg and R Hya. 
The radial velocity of the absorption lines in x Cyg 
is constant, as discussed in Sect. 5.2. The Mira S star 
S UMa (=HD 110813) is perhaps a binary, since the radial- 
velocity period P — 592.2 d (Table 8) is well distinct from 
the 225.9 d period of the Mira variations. More measure- 
ments arc needed, however, before that orbital solution 
may definitely be accepted. 

As far as Tc-poor carbon stars are concerned, X Cnc 
(= HD 76221) is probably binary, with a 491 d orbit (quite 
distinct from the -uncertain- 195 d period quoted by the 
GCVS for the SRb variations). 

A possible orbital solution, with P = 544.2 ± 5.7 d and 
e = 0.55, has been found for the SC star BD -08° 1900 
(Table 8 and Fig. 6). A light curve for this star is pro- 
vided by Jorissen et al. (1997), who find photometric vari- 
ations on a time scale of about 59 d. The absence of coin- 
cidence between the radial- velocity and photometric (qua- 
si-)periods may lend some credit to the orbital nature of 
the radial- velocity variations observed for BD —08° 1900. 
We do not accept this interpretation, however, without 
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the following word of caution: another SC star, GP Ori 
(=HD 286340), exhibits radial- velocity variations of a na- 
ture very similar to those of BD — 08°1900 (Fig. 6; note 
especially the drift observed in the early phase of the mon- 
itoring, reminiscent of that observed for BD — 08°1900), 
but no orbital solution could be obtained for GP Ori. At 
this point, we cannot exclude the possibility that the 'or- 
bital' solution found for BD —08° 1900 is simply a con- 
sequence of a favourable time sampling of irregular, in- 
trinsic variations. One should note in that respect that 
the cry/^vr ratio observed for BD -08° 1900 fits well the 
predictions of a simple linear model of adiabatic acous- 
tic oscillations (Jorissen et al. 1997). Finally, the orbital 
parameters of BD —08° 1900 would locate that star in a 
unusual region of the (e, logP) diagram (Fig. 4 of Joris- 
sen et al. 1998), which is another argument against that 
orbital solution. 

In fact, relatively large- amplitude radial- velocity varia- 
tions (<7vr of a few km s _1 ; Table 3e of Jorissen et al. 1998) 
are a common feature among SC stars. For the two Mira 
SC stars in our sample (R CMi and RZ Peg), periodic 
radial- velocity variations were found with the same period 
as the light curve (Sect. 5.1). For the remaining SC stars, 
which are of semiregular (SR) or irregular (L) variability 
types, very irregular radial-velocity variations are indeed 
observed, with the possible exception of BD —08° 1900 dis- 
cussed above. When the number of measurements is rela- 
tively small (< 10) and covers a limited time span, these 
variations may mimick orbital variations. 
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